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Abstract. - We study the kinetics of domain growth in ferromagnets with random exchange 
interactions. We present detailed Monte Carlo results for the nonconserved random-bond Ising 
model, which are consistent with power-law growth with a variable exponent. These results are 
interpreted in the context of disorder barriers with a logarithmic dependence on the domain 
size. Further, we clarify the implications of logarithmic barriers for both nonconserved and 
conserved domain growth. 



A homogeneous binary mixture becomes thermodynamically unstable if it is rapidly quenched 
below the coexistence curve. The subsequent far-from-equilibrium evolution of the system is 
characterized by the emergence and growth of domains enriched in the new equilibrium phases. 
The domain morphology is quantified by (a) the time-dependence of the domain scale R(t), 
where t is the time after the quench; and (b) the correlation function or its Fourier transform, 
the structure factor [1] . There is a good understanding of domain-growth kinetics in pure and 
isotropic systems, where the domain scale shows a power-law behavior, R(t) ~ t e . For the 
case with nonconserved order parameter, e.g., ordering of a ferromagnet into up and down 
phases, we have — 1/2. On the other hand, for the case with conserved order parameter, 
e.g., phase separation of a binary (AB) mixture into A- and B-rich domains, we have 8—1/3 
when growth is driven by diffusion. 

Of course, real experimental systems are neither pure nor isotropic. In this letter, we focus 
on domain growth in ferromagnets and binary alloys with quenched disorder in the form of 
random exchange interactions. There have been many experimental [2-4] and numerical [5-10] 
studies of this problem [11]. At early times, domain coarsening is not affected by disorder. 
Then, there is a crossover to a disorder-affected regime, which occurs earlier for higher disorder 
amplitudes. There have also been many studies of domain growth in spin glasses [12,13], where 
the amplitude of disorder is such that the local exchange coupling may be either ferromagnetic 
or antiferromagnetic. Inspite of this attention, the nature of asymptotic domain growth in 
both random magnets and spin glasses remains the subject of much controversy. The present 
letter resolves this controversy in the context of random magnets. We present detailed Monte 
Carlo (MC) results which show that asymptotic growth in these systems is consistent with 
a power-law behavior, with the growth exponent (0) depending on the temperature (T) and 
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disorder amplitude (e). Further, we present analytical arguments which clarify the functional 
form of 0(T, e) for systems with both nonconserved and conserved order parameters. 
We study the random-bond Ising model (RBIM) with the following Hamiltonian: 

H = —^^JijSiSj, Si = ±1, (1) 

m 

where the exchange couplings Jij are drawn from a probability distribution, e.g., Gaussian, 
uniform, bimodal, etc. We focus on the nearest-neighbor case, denoted by the subscript (ij). 
A kinetic version of the RBIM is obtained by associating Glauber spin-flip kinetics (Si — ► — Si) 
or Kawasaki spin-exchange kinetics (Si <-> Sj) with the Hamiltonian in Eq. In this letter, 
we consider the ferromagnetic case, where Jy > always. The case where Jij can be both 
> (ferromagnetic) and < (antiferromagnetic) is relevant to spin glasses. 

An important study of the nonconserved RBIM is due to Huse and Henley (HH) [15]. 
HH argued that coarsening domains are trapped by energy barriers E B (R) ~ E Q R^\ with 
exponent ip — x/(2 — £), where x an d C are the pinning and roughening exponents. For d = 2, 
these exponents are known to be x = 1/3 and £ = 2/3 [16], yielding ip — 1/4. For d = 3, a 
perturbative calculation gives -0 — 0.55 [15]. To obtain the corresponding growth law, it is 
convenient to use the framework of Lai et al. (LMV) [17]. For curvature-driven growth in 
nonconserved systems, the length scale obeys 

R = a(R,T)/R, (2) 

where a(R, T) is the diffusion constant which depends (in general) upon both the length 
scale and the temperature. In the presence of energy barriers, domain growth proceeds 
via thermally-activated barrier hopping with a(R,T) ~ ao exp(— Eg/T). Replacing the HH 
energy-barrier scaling in Eq. @, we obtain the crossover behavior: 

R(t)aRo(T,e)h(j^, with (3) 

' rp \ l/l/> ^ / rp \ 2/V' 



MT,e) = ( — I , *o(T,e) = — I - ) , and (4) 
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h(x) — I —xj , x <C 1, 

= (In a;) 17 "', x > 1. (5) 

The logarithmic growth law proposed by HH has motivated many experimental and nu- 
merical studies of domain growth in random magnets. However, to date, there is no clear 
confirmation of this growth law. For example, the experiments of Ikeda et al. [2] found that 
there is no universal logarithmic law. Numerical studies have also found no evidence for a 
HH-type law, valid over extended time-windows and parameter regimes. On the contrary, 
some experiments [4] and simulations [6] suggest power-law growth with a variable exponent. 

Let us next present results from our MC simulations, which were done on L x L lattices 
(in d = 2) with periodic boundary conditions. The results presented here were obtained for 
uniform -distributions on the interval [1 — e/2, 1 + e/2] or [1 — e', 1], where e and e' quantify 
the degree of disorder. The limits e = (or e' = 0) and e = 2 (or e' = 1) correspond to 
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Fig. 1 - Domain growth in the RBIM with Glauber kinetics. We show evolution pictures at t = 10 5 
MCS for a 512 x 512 lattice, after a quench from T = oo to T = 0.5. The up spins are marked in 
black, and the down spins are unmarked. The snapshots correspond to different disorder amplitudes, 
e = l,2. 

zero and maximum disorder, respectively. Assigning random initial orientations to each spin, 
we rapidly quench to T < T c ~ 2.269. All our statistical data is obtained with lattice sizes 
L = 1024, as an average over 50 independent initial conditions and disorder configurations. 
We also produced data for L = 256 and 512 (not shown) to ensure that our results are not 
influenced by finite size effects. 

Figure ^ shows typical evolution snapshots (at t = 10 5 MCS) for different disorder am- 
plitudes e = 1,2. As expected, the length scale at a given time decreases with increase in e. 
However, the domain morphology does not visually differ from that for the pure case. We 
have confirmed (not shown here) that the scaled correlation function [C(r,t) vs. r/R(t)] is 
independent of the disorder amplitude, and is numerically equivalent to that for the pure case 
(e = 0). This observation has earlier been made by Puri et al. [8] and Bray and Humayun [9]. 

Next, we consider the time-dependence of the domain size, which is obtained as the distance 
over which the correlation function decays to half its maximum value. In Fig. [21 we undertake a 
direct test of the FTFT law by first plotting, i? 1 / 4 vs. t on a log-linear scale for different quench 
temperatures and e' = 1, with Jy E [1 — e',1]. The data does not satisfy the asymptotic 
HH growth law, which corresponds to a straight line on this plot. More generally, we have 
attempted to fit the data to the logarithmic function, hit = aR x + b. This function does 
not give a good fit to the data. For acceptable fits, the exponent x depends strongly on the 
temperature (cf. Ref. [2]), in contrast with the prediction of a universal growth law. Similar 
results are obtained if we fix T and vary e, and we do not present these here. 

We also tried to fit the data sets to the crossover scaling form described by Eqs. (3)- (5), 
the result for which is shown in the right portion of Fig. 2. We record the following points of 
disagreement with the proposed scaling: (1) The short-time behavior is not described well by 
Eq. (5), i.e., in the plot g(x) does not fit well to the scaling curve. (2) The proposed asymptotic 
behavior, i.e., the curve f(x), does not fit the scaling curve well even for the largest times. 
(3) The temperature-dependence of the crossover length Rq(T) and the crossover time to(T) 
is stronger than a power law, which is incompatible with the expected behavior in Eq. (4). 
The parameters Eq and clq occurring in Eq. (4) are expected to decrease with increasing 
temperature, and therefore Rq and to may be expected to increase faster than T 4 and T 8 , 
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Fig. 2 - Left: Plot of R 1/4 vs. t on a log-linear scale for e — 1, i.e. Jy is uniformly distributed 
between [0, 1], and different quench temperatures. We show data for 5 temperature values: T = 
0.1,0.2,0.3,0.5,0.7. Right: Scaling plot according to Eqs. (3)-(5) (Jij as in left figure). For each 
temperature T the values for Ro(T) and to{T) have been chosen so as to obtain a smooth scaling curve 
h(x). The functions g(x) oc x 1/l2 and f(x) oc (lna;) 4 represent the expected asymptotic behaviors for 
i < 1 and x 3> 1, respectively, according to Eq. (5). The inset shows the temperature-dependence 
of the fit values Ro(T) and to(T) and their expected T-dependence according to Eq. (4), which is T 4 
and T 8 , respectively, if ao and Eq are only weakly (sublinearly) T-dependent. 



respectively. However, their putative T-dependence turns out to be much too strong. Note 
that the crossover time to in the inset of Fig. 2 varies over 20 decades when T varies over only 
half a decade from 0.2 to 0.7, and we do not see why the pinning energy amplitude Eq or the 
domain wall tension should have such a strong T-dependence. Based on observations 1-3, we 
believe that the crossover scaling form in Eqs. (3)- (5) does not describe the data well and we 
suggest here an alternative picture. 

Several groups [4, 6] have reported that random magnets exhibit power-law growth with 
variable exponents. In Fig. [31 (left), we plot R vs. t on a log- log scale for the data in Fig- El 
For low T, the plots are linear over 5 decades in time. For high T, the log-log plots exhibit 
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Fig. 3 - Left: Log-log plot of 7? vs. t for the length-scale data shown in Fig. Right: Log-log plot 
of R vs. t for T = 0.5 and different disorder amplitudes: e = (pure case), and e = 0.1, 0.5, 1, 1.5, 2. 
Here Jij is uniformly distributed over [1 — e/2, 1 + e/2]. 
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a crossover behavior. In Fig. [3] (right), we present a similar plot for T — 0.5, and different 
values of the disorder. Again, the data exhibits a power-law behavior for large values of e, 
and crossover behavior for small values of e. Our data is consistent with power-law growth 
with a variable exponent, at least for low T and high e. We have obtained similar results for 
domain growth in the dilute Ising model (DIM), and will present these elsewhere [21]. 

Let us understand the origin of growth exponents which depend on T and e. For the 
DIM, Henley [18] and Rammal and Benoit [19] have argued that the fractal nature of domain 
boundaries results in a logarithmic i?-dependence of trapping barriers. We propose that this 
is generally applicable [14] and examine the implications thereof. At early times and small 
length scales, we expect disorder-free domain growth. This suggests the barrier-scaling form: 

E B (R)~eln(l + R), (6) 

where R is measured in dimensionless units. Substituting a(R,T) ~ clq exp(— Eg/T) in 
Eq. J2J), we obtain 

R = a (l + R)- e / T /R. (7) 

The solution of Eq. @ is the crossover function 

R(t) ~ (2a t) 1/2 1 t«to, 

~ [(2 + e/T)a tf T ' e) , t » t , (8) 

with the growth exponent 

9(T,e) = (2 + e/T)-\ (9) 

The early- and late-time behaviors in Eq. JSJ arise in the limits i{<l and R 3> 1, respectively. 
The crossover length and time can be identified by reformulating Eq. (JSJ) as Eq. © with 

R o(T,e) = p e y} (1 _26)> *o( T > e ) = ^( 2(9 20)i/(i-20)' and ( 10 ) 



h(x) = x 1/2 , x < 1, 

= x e , x > 1. (11) 

In Fig. 21 we plot 1/6 vs. 1/T and e for the data shown in Fig. [3] The resultant linear 
plots confirm the logarithmic barrier-scaling scenario proposed above. For high temperatures 
and weak disorder amplitudes, recall that the length-scale data exhibits a crossover behavior. 
In these cases, we estimate the asymptotic exponent by tracking the effective exponent 8 c ff = 
d(lnR)/d(\nt) as a function of time. 

We note the following points: (1) When measuring the growth by alternative methods, 
we get the same estimates for the growth exponents. For instance, we applied the following 
Monte Carlo renormalization group approach. First, we scaled the lattice twice by taking 4 
block spins to reduce thermal noise. As a result, a 1024 2 lattice becomes a 256 2 lattice. Then 
we scanned horizontally and vertically through the lattice and measured the total number 
of domain walls. Finally, we divided the system area by this number to get the average 
domain length. From the time-dependence, we extracted the exponents shown in Fig. 0] 
(left), and we observe that they agree well within numerical accuracy. (2) We checked that 
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Fig. 4 - Left: Exponent 1/9 vs. 1/T for the data shown in the left part of Fig. [3] (circles), and 
obtained for the same systems using estimates for the domain size via the lattice scaling method 
described in the text (squares). Note that for T > 1 » T c (e' = 1), i.e., for 1/T < 1 the system is 
in the paramagnetic phase. Hence R cannot grow algebraically any more and no data for 1/9 are 
available in this region. Right: Exponent 1/9 vs. e for the data shown in the right part of Fig.|3] 



the correlation functions that we calculated obey the expected dynamical scaling behavior: 
C(r,t) =c[r/R(t)] ^ c[r/t 9 ^}. 

Given our results in the nonconserved case, we expect that the length-scale data for the 
conserved case should also exhibit power-law behavior with a variable exponent. We have 
undertaken MC simulations of the RBIM with conserved kinetics, and the results will be pre- 
sented elsewhere [21]. Here, we confine ourselves to discussing the implications of logarithmic 
energy barriers for growth exponents in conserved systems. 

In the absence of disorder, the domain scale obeys the Huse equation [20] R = Do/R 2 , 
with the solution R(t) ~ (3-Doi) 1 ^ 3 , which is referred to as the Lifshitz-Slyozov growth law. 
As before, the presence of disorder renormalizes the diffusion constant Do by an Arrhenius 
factor. For logarithmic barriers as in Eq. ©, the corresponding kinetic equation is 

R = D (l + R)- e/T /R 2 . (12) 

The short- and long-time solutions of Eq. H12|) are obtained as follows: R(t) ~ (SDot) 1 ^ 3 for 
t < t and R(t) ~ [(3 + e/T)D t} e{Tte) for t > t , where 

0(T,e) = (3 + e/T)-\ (13) 

The asymptotic exponent differs from that for the nonconserved case. This should be 
contrasted with the HH scenario, where the asymptotic growth law is the same for the 
nonconserved and conserved cases [1,8], which is seen by incorporating the HH barrier- 
scaling form in R — Do/R 2 . The crossover form of the solution of Eq. I|12|) is Eq. @ with 
Ro(T,e) = 1/{W)^ (1 ^ and t (T,e) = D' 1 / '(3030)1/(1-30) and h ^ = x i/a for x < 1 and 
h(x) = x for x ^> 1. 

In conclusion, we have presented detailed MC results for domain growth in random magnets 
for a wide range of temperatures and disorder amplitudes. Our results do not support the 
HH scenario of logarithmic domain growth. (Of course, we cannot rule out the possibility of 
a logarithmic regime at even later times than those investigated here.) Rather, they are in 
agreement with results [4,6] which exhibit power-law growth with a variable exponent. This 
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scenario arises naturally in the context of logarithmic energy barriers, and the corresponding 
functional dependence of 9 on T and e is in excellent agreement with our numerical results. 
The results in this letter provide a framework for the analysis of experiments and simulations 
on domain growth in disordered magnets and binary mixtures. 

* * * 
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